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Abstract 

We study the chiral symmetry structure in a linear sigma model with 

fermions in the presence of an external, uniform magnetic field in the 'ef- 
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fective potential' approach at the one loop level. We also study the chiral 
phase transition as a function of density in the core of magnetized neutron 
stars. 
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In the absence of quark mass matrix, QCD is invariant under chiral 
transformation at the lagrangian level. However, the dynamics of 
QCD are expected to be such that chiral symmetry is spontaneously 
broken with the vacuum state having a nonzero quark- antiquark con- 
densate (0 \qq\ 0} and the Goldstone theorem then requires the exis- 
tence of approximately massless pseudo-scalar mesons in the hadron 
spectrum. At high temperatures and/or at high densities, the quark 
condensates are expected to melt at some critical point and chiral 
symmetry is restored |T|. Chiral phase transition and phenomeno- 
logical consequences in the form of experimentally observable sig- 
natures have been extensively discussed in the literature j2j . It has 
also been suggested [f3], |j that systems with spontaneously broken 
symmetries may make a transition from broken symmetric to re- 
stored symmetric phase in the presence of external fields. It has 
been shown [|5], U that there exist some realistic models where the 
symmetry restoration takes place. In QED uniform, external static 
magnetic field is known to break chiral symmetry dynamically at 
weak gauge couplings |7j . 

Large magnetic fields with strengths upto 10 18 gauss have been 
conceived to exist |8j at the time of supernova collapse inside neu- 
tron stars and in other astrophysical compact objects and in the 
early Universe. Effect of such a strong magnetic field on chiral phase 
transition is thus of great interest for baryon free quark matter in 



the early universe and for high density baryon matter in the core of 
neutron stars. To study chiral phase transition in QCD we need a 
nonperturbative treatment. Lattice techniques and the Schwinger- 
Dyson equations provide specially powerful methods to study the 
chiral structure of QCD [fTQ] . Klevansky and Lemmer [|TTj stud- 
ied the chiral symmetry behaviour of hadronic system described by 
Nambu- Jona-Lasinio model minimally coupled to a constant electro- 
magnetic field. Solving the gap equation they found that whereas a 
constant electric field restores chiral symmetry, a constant magnetic 
field inhibits the phase transition by stabilizing the chirally broken 
vacuum state. This conclusion was confirmed by later studies of 



the NJL model |I2fl . Shushpanov and Smilga [jT^ considered QCD 



with the massless flavors in the leading order of chiral perturbation 
theory and studied the dependence of quark condensate on an exter- 
nal magnetic field by studying the schwinger- Dyson equation and 
showed that an external magnetic field increases the condensate. 

A particularly attractive frame work to study such systems is 
the linear sigma model originally proposed as a model for strong 
nuclear interactions [|P| . We will consider this as an effective model 
for low energy phase of QCD and will examine the chiral symmetry 
properties at finite density and in the presence of external magnetic 
field. To fix ideas we consider a two flavor SU(2) x SU(2) chiral 



quark model given by the lagrangian. 

where i/j is the quark field a and ir are the set of four scalar fields 
and g is the quark meson coupling constant. The potential U(cr, 7?) 
is given by 

U(*,*) = -±f( a ' + iP) + \\( a > + iPf (2) 

For jj? > chiral symmetry is spontaneously broken. The a field 
can be used to represent the quark condensate, the order parameter 
for chiral phase transition and the pions are the Goldstone bosons. 
At the tree level the sigma , pion and the quark masses are given by 

ml = 3AcrJ - fi 2 ; m 2 n = \o 2 cl - /i 2 ; raj = ga d (3) 

where o 2 d — \- = f 2 . An elegant and efficient way to study symme- 
try properties of the vacuum at finite temperature, density and in 
the presence of external fields is through the "Effective Potential" 
approach discussed extensively in the literature |}I5fl. We will com- 



pute here, in the one loop approximation, the effective potential in 
the presence of external magnetic field which is defined through an 
effective action V(a, B) which is the generating functional of the one 
particle irreducible graphs. The effective potential is then given by 

V eff (a,B) = V Q (a) + V 1 (a,B) (4) 



where Vi(a, B) is obtained from the propagator function G(a, B) by 
the usual relation Vi(a,B) = — ^.Tr log G (a, B). 

Alternatively one can compute the shift in the vacuum energy 
density due to zero-point oscillations of the fields considered as an 
ensemble of harmonic oscillators [|T5j . We thus require energy eigen- 



values (excitations) of particles in the magnetic field, which can be 
easily obtained, and in the absence of anomalous magnetic moment 
for uniform static magnetic field in the z-direction for a particle of 
mass M, charge q and spin J, are given by [TI| 



E{k z: n, J z ) = {k 2 + M 2 + (2n + 1 - sign(q) j z ) \q\ B) 2 (5) 

where n represents the landau level. Contribution of scalar particles 
of mass M to V\{M 2 ) after wick rotation is thus given by 

V l (M*) = 1 -J^- 4 ln(kl + M*- t e) (6) 

In the presence of magnetic field, all we need to do is to replace 
the phase space integral J 4% by || E^Lq A^p an d the energy by 
expression (0) for charged particles. For a scalar field of charge ±e, 
we thus have 

Vi(M 2 ,B) = f-tj^%ln(k 2 e + (2n+l)eB + M 2 ) 

eB d T(a-d/2) 

~ 4^ da r(a)(47r) rf / 2 
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eB ,. d T(a-I) 
lim — — 



4tt a^o da r(a)(47r) a -i 
1 ,, d M 2 + eB 



(2eB) a ~i K 2' 2eS 



a =0,(i=2 



(7) 



where ((z, q) is the generalized Riemann zeta function 

00 1 A roo fZ-l p -qb 

«*'*) = £ r^ = frr/o d * r^7 W 

«=o (9 + n) z T(z) J o 1-e t 

The potential (0) has poles at a=0, 1 and 2 for d=2 which can be 
absorbed in the counter terms. The finite part depends on the exact 
renormalization conditions that are imposed. In what follows we 



would use the MS renormalization scheme. From eqns. (0) and fl8|) 
we can write 



l-a „-& 



eB ,. 5 (2eB) L ~ a r ,. .„_, e" 



Fi(M", 5) - -^r lim ^ v ~r; y N / eft T" 2 — -r (9) 
327r 2 «-o5a r(a) •/ sinhf 



which converges for Re a > 2. We analytically continue the result in 
the complex a-plane and use dimensional regularization technique to 
extract the finite contribution. To proceed further we first consider 
the case ^ < 1, expand e _ 2Ss* and formally integrate (0) to obtain 

m 



\n\ R r) °° M 2 (— 1 Y 

v x {m 2 ,b) = -™£ lim ir £ ' w+ry_1 



32tt 2 «^o da n = 2e5 y */! (M 2 )"- 1 

2 a+v ~ 1 - 1) T(a + v - 1) C(a + y - i)(io) 
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Keeping leading terms in ^b we obtain 

Vi{M\B) = --±j[fg! C (2)log2ei3 + 

€ -^- log2-M 4 | log2eB + ..] (11) 

The leading term for the contribution of charged Goldstone bosons 
relevant for symmetry considerations is 

F 1 (M 2 , J B)^-^^log2 (12) 

which agrees with the earlier results |5|, [T3f upto a factor of order 
one. For the case of ^ > 1 we write (0) as 

Vl (M\B) = - -L, lim #- J- rdx^e-"' ^ 



32-7T 2 «-»-o 5a r(a) -'o sinhe^Bx 

(13) 

and keeping leading terms obtain 

^(M 2 ,^) ~ ^[M 4 (logM 2 - 5) - ^(e5) 2 logM 2 ] (14) 



which agrees with the result obtained by Salam and Strathdee 
Likewise for the charged fermion fields using (|5]) we obtain 

Vi M 2 ,5 = ^Vlim — v ^ / dt ^"V^coth- 

v ' ; 32tt 2 a^o 5a r(a) ^o 2 

(15) 

The factor of 4 and positive sign account for the spinor nature of 

the fermi field. In the limits mentioned above, we obtain 

\n\ RM 2 
F X (M 2 , B) = 1 \ f 2-(l " logM 2 ) (16) 
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and 

^(M 2 , B) ~ "^[M 4 (logM 2 ~\) + \ (\q\ BflogM 2 ] (17) 

for -fjr- > 1 and < 1 respectively. 

The total V e ff(a, B) for the sigma model at the one loop level is 
thus given by 



2 ,,2 



1 /ox 2 2x2 , ,3A^- M ' 3, 



(3Af7 2 - ^ ) log 



647r 2 ' m 2 2 
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(/;t (Aa 2 -/i 2 ) log 2 



16tt 2 



N c 44 g 2 a 2 3, 

ElflM lo g 



!6tt 2 /i, L " " v "'"' m 2 2 

+ 1 d^i^) 2 iog C5" ] (i8) 

For ^p > 1 ,the last term in eqn.( |I8|) is replaced by eq. fllfi| ) 
summed over flavors. In figure 1, we plot V e ff(a,B) as a function 
of a for different values of magnetic field and compare it with the 
case of zero magnetic field by ignoring the B independent one loop 
terms. As input parameters we choose the constituent quark mass 
ruf— 500 MeV, sigma mass m a — 1.2 GeV and 7^=93 MeV. We find 
that in the presence of intense magnetic fields the chiral symmetry 



breaking is enhanced. For magnetic field large compared to m^ , 
from eqn. ( p70|) we observe that though the fermionic contribution is 
towards symmetry restoration, it is not enough to offset the contri- 
bution of charged goldstone pions. 

In order to study chiral symmetry restoration in the case of neu- 
tron stars as a function of chemical potential fi associated with finite 
baryon number density we employ the imaginary time formalism by 
summing over Matsubara frequencies. This amounts to adding the 
fermionic free energy to the one loop effective potential and is given 

by 

^) = -J/(^ln(l + e-«*-")) (19) 

which in the presence of static uniform magnetic field becomes 

where 7 is the degeneracy factor and is equal to 2N C for each quark 
flavor. We consider cold dense isospin symmetric quark matter for 
which the integrals can be performed analytically. The baryon num- 
ber density corresponding to the chemical potential \x is given by the 
usual thermodynamical relations. 

WMl^E ^(^-^V) 1 (21) 

6 flav O 71 " 



and 



1 71 II £? 

Nb{^ B) = - £ 1J ty~ ( 2 " <W) VV " 3 2 * 2 ~ 2n M B (22) 
6 „. = n 47T Z 



n=0 
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for zero and finite magnetic field respectively. Here n max = Int 
To study chiral symmetry behavior at finite density in the presence 
of uniform magnetic field, we minimize effective potential with re- 
spect to the order parameter a for fixed values of chemical potential 
and magnetic field ( which then fixes the baryon density ). The 
results are shown in figure 2 where we have plotted the order pa- 
rameter a as a function of density at T=0 for different values of 
magnetic field. The solution indicates a first order phase transition. 
The actual transition takes place at the point where the two minima 
of the effective potential at cr=0 and a=a(/i,B) non zero become 
degenerate. The lower values of o (shown by dotted curves) are 
unphysical in the sense that they do not correspond to the lowest 
state of energy. We find that magnetic field continues to enhance 
chiral symmetry breaking at low densities as expected but as the 
magnetic field is raised the chiral symmetry is restored at a much 
lower density compared to the free field finite density case. This 
can be clearly seen from figure 3 where we have plotted the phase 
diagram in terms of baryon density and magnetic field. 

In conclusion we have examined the chiral symmetry behavior of 
the Linear Sigma model in the presence of static, uniform magnetic 
field at the one loop level at zero density and at densities relevant 
in the core of neutron stars. We find that the contribution of scalar 
and fermion loops leads to an increase in chiral symmetry breaking. 
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At high densities too, this effect persists and for magnetic fields of 
strength upto 10 18 Gauss, there is enhancement in chiral symmetry 
breaking resulting in the restoration of symmetry at densities higher 
than if no magnetic field were present. However, in the case of high 
magnetic field B > 10 19 Gauss the chiral symmetry is restored at 
lower densities. Thus in the core of neutron stars, if the nuclear mat- 
ter undergoes a transition to deconfined quark matter, the presence 
of magnetic field would imply the existence of massive quark matter 
due to enhancement in chiral symmetry breaking. This would affect 
the equation of state and will have astrophysical implications. 

Acknowledgements 

We would like to thank Professor J.V. Narlikar for providing hos- 
pitality at the Inter-University Centre for Astronomy and Astro- 
physics, Pune 411 007, India where this work was initiated. 



11 



Figure captions 

Figure 1. Effective Potential in units of (lOOMeF) 4 as a function 
of a(MeV) for different values of the magnetic field. The curves 
a, b, c, d and e are for B=0, 10 16 , 10 18 , 10 19 and 3 x 10 19 Gauss 
respectively. 

Figure 2a. Chiral condensate a(MeV) as a function of baryon 
density in f m 3 for magnetic field B=0. 

Figure 2b. Chiral condensate a(MeV) as a function of baryon 
density in f m 3 for magnetic field B=10 16 Gauss. 

Figure 2c. Chiral condensate a(MeV) as a function of baryon 
density in f m 3 for magnetic field B=10 18 Gauss. 

Figure 2d. Chiral condensate a(MeV) as a function of baryon 
density in f m 3 for magnetic field B=10 19 Gauss. 

Figure 3. Phase diagram as a plot of magnetic field versus baryon 
density in units of frrT d . 
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